Capillary-gravity wave resistance in ordinary and magnetic fluids 
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Abstract 



i 

Q ', Wave resistance is the drag force associated to the emission of waves by a 

I moving disturbance at a fluid free surface. In the case of capillary-gravity 

waves it undergoes a transition from zero to a finite value as the speed of 
Q_i' the disturbance is increased. For the first time an experiment is designed 

in order to obtain the wave resistance as a function of speed. The effect of 
viscosity is explored, and a magnetic fluid is used to extend the available 
1^ I range of critical speeds. The threshold values are in good agreement with the 

fSJ ' proposed theory. Contrary to the theoretical model, however, the measured 

■ wave resistance reveals a non monotonic speed dependence after the threshold. 

O ; PACS numbers: 47.35. +i, 68.10.-m, 75.50.Mm 
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When an object is moved at the free surface of a fluid, it experiences a drag force 
which physically originates from: (a) bulk dissipation in a viscous boundary layer for low 
Reynolds numbers, and in the turbulent wake for high Reynolds numbers; (b) the emission 
of capillary-gravity surface waves. Such waves remove momentum from the perturbating 
object to infinity. The associated force that the object experiences is called wave resistance. 
For the convenient moderate speeds on which we focus in this paper it may overcome the 
bulk dissipation type drag. 

Wave resistance has been studied for more than a century in the case of pure gravity 
waves |]T|, mainly because this topic has obvious naval applications 0]. In this case the 
variation of the wave resistance R with the speed V follows two regimes. A critical velocity 



ygrav jg jj^iposed by the characteristic size L of the ship: V^'^'^^ = ^gL/2n, g being the 
gravity acceleration. For V < V^^"''" the wave resistance is very close to zero and behaves as 
R oc W - Vf'"-" for V > 1//™''. This has recently been analyzed in terms of an imperfect 
bifurcation by one of the co-authors 0. 

The case of capillary-gravity waves has been theoretically treated in a recent work of 
Raphael and de Gennes @]. Such waves are generated when the size of the perturbating 
object is small compared to the capillary wavelength Ac = 27i{a/ pgY^^, where a is the surface 
tension of the free air-fiuid interface and p the density of the fluid. The dispersive properties 
of capillary-gravity waves are such as there exists a minimum phase speed Vc = {Aag/ pY^^ 
at which waves are able to propagate. Since the pattern is stationary in the reference frame 
of the moving object, no wave can be emitted for V < \4 [Q, and therefore there is no wave 
resistance in that case. As it has been shown in [Q, the wave resistance experiences a finite 
jump Rc aX V = Vc and increases above Vc- The system is thus supposed to undergo a 
discontinuous bifurcation. 

In order to check these theoretical predictions it is necessary to vary Vc by means of p 
and (J variations. By adding a surfactant to water, a may be easily chosen between say 
20 and 73 mN/m. Consequently Vc will merely vary from 17 cm/s to 23 cm/s. A more 
efficient control parameter is thus needed. We here show that the action of a magnetic 
field on a magnetic fiuid provides a means to tune the critical velocity from its maximum 
value V^^^ = {Aag/py/^ down to 0. Using a magnetic fiuid, along with other regular 
fiuids of different viscosities, we perform RiV) measurements, as the problem has not been 
experimentally treated yet. 

In a regular fluid, the wave emission process is controlled by the dispersion equation of 
capillary-gravity surface waves, u'^ = gk + ak^/p, where u is the circular frequency and k 
the modulus of the wave vector. The condition for stationarity of the wave pattern in the 
frame of reference of the moving object is = kVcos6 , where 6 is the angle between the 



speed and wave vectors. Thus we obtain the following equation {kc = \J pg / a is the capillary 
wavevector): 






(1) 



which has no solution for V <Vc- For a moving Dirac Delta pressure distribution P(x, y, t) 
pS{x — Vt,y) the wave resistance is 0]: 
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where k^{6) and k^{6) are the two roots of Eq. (|l|). This formula remains vahd as long as 
the characteristic size of the pressure distribution in experiments is much smaller than the 
capillary wavelength. In those conditions, close and above the threshold, the wave resistance 
has a finite value Rc = p'^kc/2\/2a and increases monotonically with speed (see inset of Fig. 
|1| and the uppermost curve in Fig. 

In order to measure R{V) for the various fluids, they are placed into a circular channel 
dug into a Teflon covered aluminum dish. The latter is fixed to a shaft and rotated at 
constant rate, thus simulating a steady flow for the fluid. The radius of the channel is 
20 cm, its width is 2 cm. A 4 cm wide channel is also used, showing no significant difference 
in the experimental results. The depth of the fluid is usually more than 1 cm, ensuring the 
validity of the infinite depth approximation. 

The disturbing object consists of a vertical bronze wire (radius r = 0.2 mm) whose tip 
just touches the surface of the fluid (the wire is wetted by a few tenths of millimeters of 
fluid). The deflection of the wire is proportional to the horizontal force exerted on its free 
end (which is typically in the order of a micronewton) . It is measured with an infrared 
optical sensor. The calibration of the sensor is obtained by tilting the base to which the 
wire is attached. A more detailed description of the measuring method will be published 
later. 

Though no theory includes 3D viscous effects so far, we measure the wave resistance 
for different viscosities. To this purpose several mixtures of water and glycerol are used: 
the surface tension a and the densities p of the mixtures are very close to one another (see 
Table 1) so that the impact of viscosity alone may be monitored in our experiments. The 
viscosities are measured with a standard Poiseuille viscometer. 

Fig.0 displays the variation of the experimental drag R^xp as a function of speed for the 
various mixtures. All the measurements are obtained by increasing and then decreasing the 
speed: there is no hysteresis. We may note that: 

(a) There is a critical velocity at which the measured drag is discontinuous. That point 
validates an important feature of the Raphael and de Gennes' theory [Q. Besides, it has 
been checked that the drag discontinuity occurs at the very speed at which the wave pattern 
develops. The measured critical velocity is 23 ± 0.5 cm/s for pure water. It corresponds 
to a surface tension interval of [65.1;77.7] mN/m into which lies the tabulated value of 
pure water surface tension 72.75 mN/m at 20 °C. For water /glycerol mixtures we obtained 
Vc ~ 22.5 cm/s, that is compatible within experimental error bars with the surface tension 
of the mixtures (around 70 mN/m). 

(b) The experimental drag is not null below the critical velocity, all the more since the 
viscosity is high. The viscous drag Rdrag that is exerted over the immersed wire must be 
added to the wave resistance R to account for the measured drag Rexp- At a speed of 10 
cm/s, the Reynolds number based on the approximate length of the wetted part of the wire 
h ~ 0.4 mm is already equal to 40 (for water). A crude estimate of the viscous drag can 
still be given by the Stokes formula : Rdrag ~ GnrihV. We experimentally check that the 
viscous drag is indeed proportional to h; it can be seen in Fig.|I| (dotted lines) that it is 
also proportional to V, at least for moderate enough speeds. The linear Rdrag dependence 
on viscosity r] is harder to assess because it is impossible to impose exactly the same h 
from an experiment to another, the wetting of the wire being imperfect; nevertheless it is 
linear within the uncertainties over h. Inset of Fig.|I| presents the R{V) variations after 
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subtraction of the viscous drag Rdrag for each sample. It is this quantity that has to be 
compared with the theoretical expression (full line in the inset). It is clear from the inset 
that a pretransitional effect takes place, as the measured drag sharply increases just below 
the threshold (the higher the viscosity, the stronger the effect). A recent model for 2D 
viscous wave resistance predicts such a feature. 

(c) The amplitude of the wave resistance discontinuity at Vc compares well with the theory. 
Assuming a perfect wetting of the wire by the fluid, the total force acting on the fluid is 
p = 2'Kra (r is the radius of the wire). Thus an estimate of the wave resistance at the 
threshold is given by Rc = ii'^r'^yJ2pga. A comparison between expected values and what is 
observed is given in Table |. The discrepancy is partially due to the imperfect wetting of 
the fluid on the wire, which leads to overestimate the applied vertical force. On the other 
hand the drag values close to the threshold fluctuate a lot. 

(d) The drag is a non-monotonic function of speed for V > Vc. In fact, it can be seen in 
inset of Fig. |l] that for > the wave resistance R first decreases as the speed increases, 
and then increases again for high enough speeds. This feature is not predicted by the current 
theory, which anyway overestimates the actual drag. Moreover it is unlikely a viscosity effect 
since it is as much marked as the viscosity is low. It is possibly a general feature of such 
a flow, and in this case the theory should be revised to include viscosity and non-linear 
aspects. 

In a magnetic fluid the dispersion equation of capillary-gravity surface waves is modifled 
with allowance for a vertical uniform magnetic fleld ||11|| : 



u; =gk + fio— — , (3) 

p Pr{Pr + l) P 

where H is the magnetic fleld, pr the relative magnetic permeability of the magnetic fluid 
(assumed to be constant |11[]) and po the vacuum magnetic permeability. For a given wave 
vector, an increase of the fleld intensity lowers the frequency of the waves. The frequency 
drops to zero when H reaches a certain critical value if* deflned by: 

[p - 1)2 Pq 

For H > the surface becomes unstable: the Rosensweig instability, sometimes called 
the peak instability, develops yielding an hexagonal array of peaks |^. The condition for 
stationarity implies that k must be a solution of: 

r X 2\ / L \ 

= 0. (5) 
Real solutions exists if and only if: 

V > V" with Vj^ = K. 




therefore a steady vertical magnetic fleld should allow the tuning of the critical velocity at 
which waves (and wave resistance) appear. 
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The wave resistance, following Eq.(0) and Eq.(|^), is given by the following integral: 



T)^ h /"arccos ^^^j— , , , , 

R»(y) = ^ / cos^ (2B{V, 6)^ + B{V, 6)-^) dd, (7) 
TTcr Jo ^ ' 

where = (^(^^ COS ^) +(;!")) 

Just above the threshold, the wave resistance has the finite value : 

Another experiment is conducted using a water based magnetic fluid synthesized accord- 
ing to the Massart method 0. Its critical field is 9.15 kA/m and its surface tension of 60 
mN/m doesn't depend on the magnetic field. Other caracteristics are given in table |I[ The 
experimental critical values and = RexpiY^) — RdragiY^) are both plotted versus 
the normalized magnetic field H/ in Fig. ^ and Fig. |^, and are compared to theoretical 
predictions @ and (H). 

The theoretical expression of (Fig- 0) remarkably fits the data points — note that 
there is no adjustable parameters. A data point lies outside the curve, but this is probably 
related to an imperfect magnetic wetting phenomenon. As the magnetic field gets closer to 
the peak instability threshold value H^,, the fluid "climbs" onto the wire, producing a much 
higher viscous drag, a situation which gets away from our inviscid linear theoretical analysis. 
This also explains the discrepancy in Fig. ^ between experimental and theoretical Rq values. 
We do not account for the force that the magnetic field is exerting at the meniscus close 
to the wire. Indeed, the very shape of the meniscus creates a non homogeneous magnetic 
field which results in a force that sucks the magnetic fluid up and changes the shape of the 
meniscus. Only advanced numerical simulations would allow to compute the net force added 

Fig. ^ presents the results obtained for different magnetic fields in a reduced representa- 
tion R/R^ = fiy/V^) with R^ = Rexp — Rdrag- It also gives a comparison to the theoretical 
predictions of Eq.(|^). As it was pointed out with regular viscous fluids, the theoretical vari- 
ations of R/Rc lie above the data points, except for H ^ H^, . Then the experimental data 
and the theory are very comparable. The present theoretical description thus gives a correct 
general trend for the influence of the field on the wave resistance. 

In conclusion, for the first time a capillary-gravity wave resistance measurement is per- 
formed on fluids of various viscosities. A drag discontinuity is always observed for a critical 
velocity Vc- Thanks to a magnetic fluid the critical velocity range is experimentally ex- 
tended. In all cases the measured critical velocities and the critical values of the resistance 
are in good accordance with the developped models. If an inviscid theory is correct at the 
threshold, there are some discrepancies for V > Vc such as a non-monotonic behavior of the 
wave resistance. Viscosity and non linear aspects should be taken into account in further 
works. Finally, in order to get rid of the viscous drag that is always present in our experi- 
ments, another mode of disturbance is to be envisaged, such as a small magnet placed just 
above the free surface of a flowing magnetic fluid. 
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FIGURES 



FIG. 1. Experimental drag Rexp as a function of speed V for different water glycerol mixtures. 
For readability purpose only a few error bars are plotted. Dotted lines : linear viscous drag Rdrag 
for each viscosity. Inset : wave resistance R = Rexp — Rdrag as a function of V . Same symbols as 
in the main figure. Full line : theoretical expression from Eq. (2). 

FIG. 2. Reduced critical speed /V^^^ at which wave resistance appears in function of the 
applied reduced magnetic field H/H^. The straight line represents the theoretical law given by 
Eq.@. There is no adjustable parameter. 

FIG. 3. Drag at threshold R^ as a function of the reduced magnetic field H/H^,. The full line 
represents the theoretical law as given by Eq.(|9|). There is no adjustable parameter. 



FIG. 4. Measured wave resistance R^ = R^xp ~ ^drag ^ function of reduced speed V/V^ 
for different reduced magnetic fields H/H^f. The theoretical curves are derived from Eq.(|^,^). The 
uppermost curve describes the wave resistance of a regular non-magnetic fluid. 
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TABLES 



TABLE I. Experimental drag discontinuity at the threshold compared to the theoretical pre- 
dictions of [4], for various water glycerol mixtures and an aqueous magnetic fluid (MF). 



Glycerol mass fraction (%) 


60 


44.5 


30 





MF 


Viscosity (mPa.s) 


12.5 


5.1 


2.6 


1.0 


7.0 


Density (g/cm^) 


1.16 


1.13 


1.09 


1.00 


1.56 


Theory (//N) 


3.9 


3.8 


3.8 


3.7 


4.2 


Experience (/xN) 


2.9 


2.6 


4.0 


3.6 


4.0 


Uiic(M;(aiiity (//X) 


0.3 


0.3 


0. 1 


l.cS 


1.0 
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